Interacting Growth Walk on a honeycomb lattice 
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The Interacting Growth Walk (IGW) is a kinetic algorithm proposed recently for generating 
long, compact, self avoiding walks. The growth process in IGW is tuned by the so called growth 
temperature T' = 1/(A:b/?'). On a square lattice and at T' = 0, IGW is attrition free and hence 
grows indefinitely. In this paper we consider IGW on a honeycomb lattice. We take contact energy, 
see text, as e = — |e| = —1. We show that IGW at /3' = cxd (T' = 0) is identical to Interacting 
Self Avoiding Walk (ISAW) at /? = ln4 (/tsT = l/ln4 = 0.7213). Also IGW at = (T' = oo) 
corresponds to ISAW at /3 = ln2 (fesT = l/ln2 = 1.4427). For other temperatures we need to 
introduce a statistical weight factor to a walk of the IGW ensemble to make correspondence with 
the ISAW ensemble. 
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We shall be concerned with a linear homo polymer 
chain, modelled by a lattice Self Avoiding Walk (SAW). 
Let z denote the coordination number of the lattice. An 
ensemble of SAW is generated by the following simple 
non-reversing blind ant algorithm. A blind ant starts at 
a site, say origin. The ant steps into one of the z nearest 
neighbour sites with a probability p — 1/z. Since the 
ant never reverses its step, in the second and subsequent 
steps, it moves to one of the z — 1 nearest neighbour sites 
with a probability (z — 1)"^. If the ant finds that the 
site has been visited earlier, then the walk is terminated 
and we start all over again. Sample loss due to violation 
of self avoidance is called the problem of attrition be- 
cause of which growing of large number of long polymer 
chains becomes difficult. Thus, a walk having N step is 
generated with a probability given by, 
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The important point is that as per the above non- 
reversing blind ant algorithm, all the N step SAWs are 
generated with the same probability. Note that the only 
interaction present is due to the excluded volume effect 
(the hard core repulsion) . In the random walk model this 
is taken care of by the self avoidance condition. 

Let us now switch on a weak interaction. The aim is 
to model the interaction that ensues when a segment of 
the polymer chain comes close to another segment. The 
segment - segment interaction can be attractive or repul- 
sive depending on the nature of the monomers present 
in the polymer chain. This interaction is usually mod- 
elled, see for e.g. [Q as follows. We say that in an SAW 
configuration, every pair of occupied nearest neighbour 
sites but not adjacent along the walk, carries an energy 
e. We call such a pair, as giving rise to a single non- 
bonded nearest neighbour (nbNN) contact. If e < then 



the interaction is attractive and if e > 0, it is repulsive. 
Thus a walk C belonging to the SAW ensemble has an 
energy E{C) — nNj\j{C) x e, where npfN{C) denotes the 
total number of nbNN contacts present in C. 

The probability with which an SAW would be found in 
a canonical ensemble of Interacting Self Avoiding Walks 
(ISAW) at temperature T = l/[kB(3] can now be ex- 
pressed as, 
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In the above, the denominator is the canonical partition 
function, given by. 



g(/3,7V)=^exp [-/3i?(C)] 
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from which all the required macroscopic properties of the 
polymers can be calculated either analytically or numer- 
ically, in the thermodynamic limit of A'^ ^ oo. In what 
follows, we consider only attractive interaction, repre- 
senting weak Van der Walls forces and set e — — |e|, see 
for e.g. [|| and without loss of generality take the strength 
of interaction as unity, i.e. \e\ = 1. 

At very high temperatures, T oo (/3 ^ 0) the nbNN 
contact interactions are unimportant; a polymer configu- 
ration is completely determined by the excluded volume 
repulsion (present always and taken care of by the self 
avoidance condition) and entropy. We get a relatively ex- 
tended configuration representing a polymer under good 
solvent conditions jsj. As T is lowered, the contact in- 
teractions become more and more important and below 
a critical temperature, called the theta point, there is 
an abrupt transition to a collapsed phase. At the theta 
point itself, we get an intermediate phase called the theta 
polymer. 
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A major problem with SAW is attrition, because of 
which we are not able to grow large number of long poly- 
mer configurations within meaningful computer times. 
Several algorithms have been proposed and investigated 
addressing the problem of attrition. These include the 
True Self Avoiding Walk(TSAW) g, Kinetic Growth 
Walk(KGW) §, Smart Kinetic Walk(SKW) @, Inter- 
acting Oriented Self Avoiding Walk (lOSAW) § etc., 
to name a few. The most recent addition to this list 
is the Interacting Growth Walk (IGW) proposed by 
Narasimhan aZ 101. 
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FIG. 1. IGW on a honecomb lattice. The first step from 
O is taken with a probability 1/3. All the subsequent steps 
until site A are taken each with a probability 1/2. If the 
walk at site A goes to site B, then an nbNN contact de- 
velops lowering the energy. The probability for this step 
is p{A ^ B) ^ exp(/3')/[l -I- exp(/3')] On the other hand 
the step from A to C does not change the energy and 
p{A ^ C) = 1/[1 + exp(/3')]. If the walk goes from A to 
B, then it is immediately followed by a step from B to D 
with unit probability. 

In IGW, we define a growth temperature T' that con- 
trols locally the growth process as described briefly be- 
low. We consider a non reversing and myopic ant rather 
than a blind ant. The probability of moving to a site 
that makes, say /i number of nbNN contacts is made 
proportional to exp(/?'/i). This is illustrated in Fig. 1, 
where we consider SAW on a honeycomb lattice. A typ- 
ical SAW trail is marked by thick connected line seg- 
ments in Fig. 1. For the honeycomb lattice z = i 
and /i = 0, 1. The walk at site marked A in Fig. 
1, can move either to the site B or to the site C. If 
the site B is selected, the energy would be reduced by 
one unit since this choice results in an nbNN contact. 
Hence we choose the site B with a relatively higher 
probability determined by the growth temperature T': 
p{A -> B) = exp(/?')/[l + exp(/3')]- On the other hand. 



the step from site A to site C, does not lead to any change 
in the energy, and p{A — > C) = 1/ [1 -I- exp(/3')]. Note 
that whenever a contact making site is selected, it is fol- 
lowed immediately by a step to the only remaining near- 
est neighbour site, and the probability for this is unity. 
In Fig. 1 , if the walk reaches site B then then in its next 
step it goes to site D with a probability p{B — > D) = 1. 
Let us say that in an N step SAW, grown as per the my- 
opic ant IGW algorithm, there are n^M contact steps ac- 
cepted and n'j^ijq contact steps avoided during the growth 
process. For example in Fig. 1,. if the step A to B is 
taken then we say the contact is accepted. On the other 
hand if the step A to C is taken then we say the contact 
is avoided. We then have. 
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1 + exp(/3' 
Let us consider two extreme cases 
Case I : 13' -> oo {T' 



For this (C) = V C, since whenever a con- 

tact making step is available, the random walk takes 
it with probability unity. In Fig. 1, this corresponds 
to p{A B) = 1 and p{A ^ C) = 0. Identifying 
E{C) = —nj\jN, we have. 
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We immediately see that the IGW ensemble at growth 
temperature T' 0' (/?' — > oo) corresponds to ISAW 
ensemble at ksT = l/ln4 = 0.7213 (/? = In 4). 

Case 2 : = (T = oo) 

At very high growth temperatures (T' — > oo), the 
nbNN contact interaction is unimportant and the walk 
steps into one of the unvisited nearest neighbour sites 
with equal probability. We have, 
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Thus at very high growth temperatures (T' oo) 
(/3' 0) IGW ensemble is equivalent to ISAW ensem- 
ble at kBT= l/ln2 = 1.4427 (/3 = ln2). 

Let us now investigate what happens when < /?' < 
oo. We rewrite Eq. (0) as. 
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The bias arising due to the contact sites avoided dur- 
ing the growth process can be removed by attaching a 
weight W{(3\ C) to a walk C generated as per the IGW 
algorithm. It is given by, 



W{P',C) 



1 + exp(/3') 



,{C) 



(9) 



The weighting for the bias removal is implemented as 
follows. Start an IGW with a weight W = 1. Ev- 
ery time a contact is available but not taken by the 
random walk, multiply the weight by a factor given by 
[1 + exp(/3')]/2. Let W{f3',C) be the weight at the end of 
an A^-step walk. Thus the IGW ensemble at (3' defined 
by the set of weights {W{P\C)} is equivalent to ISAW 
at (3 = lnF(/3'), where is given by Eq. (||). 

Infact, if we accept a weighted IGW ensemble, we can 
choose the weights appropriately (called importance sam- 
pling 1^) to make a correspondence with ISAW ensemble, 
such that (3' of IGW is the same as (3 of ISAW. This is 
carried out as follows. Let, 



exp(/?')' 



(10) 



Then we have, 



= exp [-P'EiC)] 
X [Fi(/3')]"""^'^' 
X [F2(/3')]"""^'^' X Vff^, 
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where E{C) = -nNN{C) and Vf;^"" is given by Eq. ^) 
with z = 3. The factors which depend on Fi and F2 can 
be taken care of by attaching suitable statistical weights 
to the configuration. Essentially, we start the random 
walk with a statistical weight W = 1. At any stage of 
the growth process, if both the nearest neighbour sites 
are not contact making sites, then we choose one of them 
with equal probability and proceed. We do not do any- 
thing to ly. In Fig. 1, when the walk goes from the site 
E to A, we do not adjust W. On the other hand if the 
step leads to a contact (the step A to B in Fig. 1), then 
W is multiplied by l/i^i(/3'). If the step avoids a contact 
(A to C in Fig. 1), then we multiply W by l/i^2(/5')- 
Note that these weight adjustments are carried out in 
the same spirit as the PERM B algorithm proposed by 



Grassberger We carry out the above weight multi- 
plication everytime a contact making step or a contact 
avoiding step is made during the growth process. Let 
W{l3',C) be the weight at the end the N step walk C. It 
is immediately seen that if the IGW is weighted as per 
the above procedure, then the growth temperature T' is 
the same as the temperature T of the canonical ensem- 
ble of ISAW. Thus we have an IGW ensemble defined by 
the weights {W{C,f3)}. We can employ this ensemble to 
calculate the macroscopic properties like energy, specific 
heat (fluctuations in energy) , end - to - end distance etc. 

In conclusion, we have shown that IGW at /?' = 00 is 
identically equivalent to ISAW at /3 = In 4. Also IGW at 
P' = is identically equivalent to ISAW at /? = ln2. For 
other values of /?' we introduce a weight factor to correct 
for the bias arising due to the contacts avoided during 
the growth process and show that the weighted IGW at 
/3' is equivalent to ISAW at (3 given by lnF(/3'), where 
F{[3') is given Eq. (H). Infact by choosing appropriate 
weights (importance sampling) for the accepted as well 
as avoided contacts we can render weighted IGW at (3' 
as equivalent to ISAW at (3 — (3' . 
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